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EXECUTIVE  SUMMARY 


The  influence  of  material  configuration  and  impact  parameters  on  the  damage  tolerance 
characteristics  of  sandwich  composites  comprised  of  carbon-epoxy  woven  fabric  facesheets  and 
Nomex  honeycomb  cores  was  investigated  using  empirically  based  response  surfaces.  A  series 
of  carefully  selected  tests  were  used  to  isolate  the  coupled  influence  of  various  combinations  of 
the  number  of  facesheet  plies,  core  density,  core  thickness,  impact  energy,  impactor  diameter, 
and  impact  velocity  on  the  residual  strength  degradation  due  to  normal  impact.  The  ranges  of 
selected  material  parameters  were  typical  of  those  found  in  common  aircraft  applications. 
Quadratic  response  surface  estimates  of  the  compressive  residual  strength  as  a  continuous 
function  of  material  system  and  impact  parameters  did  not  correlate  as  well  as  impact  resistance 
because  of  bifurcations  between  failure  modes.  For  a  fixed  set  of  impact  parameters,  regression 
results  suggest  that  impact  damage  development  and  residual  strength  degradation  is  highly 
sandwich  configuration  dependent.  Increasing  the  number  of  facesheet  plies  and  the  thickness  of 
the  core  material  resulted  in  the  greatest  improvement  in  the  damage  tolerance  characteristics. 
An  increase  in  the  impactor  diameter  and  impact  energy  results  in  a  significant  decrease  in  the 
estimated  residual  strength,  particularly  for  those  sandwich  panels  with  thicker  facesheets.  The 
effects  of  impact  velocity  on  damage  formation  and  loss  of  strength  were  also  addressed.  The 
development  response  surfaces  for  the  compressive  residual  strength  may  have  limited  use, 
particularly  as  the  measured  data  at  low-energy  impact  levels,  small  impactor  diameter,  and  thin 
facesheets  exhibit  significant  scatter.  Employing  the  methodology  outlined  here,  it  may  be 
possible  to  tailor  sandwich  composite  designs  to  obtain  enhanced  damage  tolerance 
characteristics  over  a  range  of  expected  impacts.  Such  efforts  may  facilitate  sandwich  panel 
design  by  establishing  relationships  between  sandwich  configuration  and  impact  parameters  that 
lead  to  improved  damage  tolerance  and  resistance. 
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1.  INTRODUCTION. 


Sandwich  construction  composites  are  used  in  a  wide  variety  of  structural  applications  largely 
because  of  their  relative  advantages  over  other  structural  materials  in  terms  of  improved  stability 
and  weight  savings.  Sandwich  plates  and  shells  are  multilayered  structures  consisting  of  one  or 
more  high-strength,  stiff  layers  (e.g.,  laminated  facings)  that  are  bonded  to  one  or  more  flexible 
layers  (core),  as  shown  in  figure  1 .  While  the  initial  design  of  composite  sandwich  structures  is 
at  a  fairly  mature  stage  of  development  [1-3],  less  progress  has  been  made  in  understanding  the 
effect  of  adverse  in-service  impact  events  on  the  structural  integrity  of  such  structures.  Such  an 
understanding  is  critical  if  there  is  to  be  widespread  use  of  sandwich  composites  in  applications 
where  structural  durability  and  damage  tolerance  are  primary  considerations  (e.g.,  aerospace  or 
automotive  applications).  Low-energy  impacts  may  induce  localized  damage  in  sandwich 
composites  (fiber  breaks,  resin  cracking,  face  sheet-core  delamination,  core  crush,  puncture,  etc.) 
and  can  be  attributable  to  a  number  of  fairly  common  discrete  sources  (hail,  tool  drops,  runway 
projectiles,  bird  strikes,  or  other  unintentional  impacts).  Any  reference  to  impact  damage  used 
herein  will  suggest  the  damage  associated  with  relatively  low-energy  and  low-velocity  impact 
events.  The  effect  of  intrinsic  processing  induced  defects  (porosity,  voids,  small  disbonds,  etc.) 
and  catastrophic  damage  associated  with  high-energy  impacts  resulting  from  airplane  crashes  or 
similar  events  are  not  considered  here. 

Characterizing  the  thermomechanical  response  of  sandwich  composites  with  varying  levels  of 
impact  damage  is  a  crucial  issue  in  the  development  of  a  damage  tolerance  plan  for  composite 
structures.  Foreign  object  impact  damage  in  sandwich  composites  may  result  in  drastic 
reductions  in  strength,  elastic  moduli,  and  durability  and  damage  tolerance  characteristics.  An 
overview  of  recent  investigations  into  the  mechanics  of  damaged  composites  may  be  found  in 
references  4-13.  The  motivation  for  this  study  is  the  development  of  empirically  based  models 
(i.e.,  “response  surfaces”)  that  isolate  the  influence  of  key  sandwich  configuration  parameters 
(e.g.,  number  of  facesheet  plies  or  facesheet  thickness,  core  thickness,  and  core  density)  and 
impact  parameters  (e.g.,  impact  energy,  spherical  impactor  diameter,  and  impact  velocity)  on  the 
damage  tolerance  characteristics  of  sandwich  composites.  This  effort  extends  the  work  of  Lacy 
et  al.  [13],  which  developed  response  surfaces  for  assessing  the  impact  damage  resistance 
properties  of  sandwich  composites.  This  study  represents  a  subset  of  a  somewhat  larger 
investigation  into  the  durability  and  damage  tolerance  properties  of  sandwich  composites  with 
varying  levels  of  impact  damage  [4,  5, 12,  and  13].  The  focus  of  this  examination  is  on  sandwich 
composites  comprised  of  plain  weave  carbon-epoxy  facesheets  and  Nomex  honeycomb  cores 
subjected  to  compression  after  impact  (CAI)  loading. 


FIGURE  1.  TYPICAL  SANDWICH  COMPOSITE  CONSTRUCTION  FORM 
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2.  DAMAGE  TOLERANCE  CHARACTERIZATION  USING  EXPERIMENTALLY 
DETERMINED  RESPONSE  SURFACES. 


In  this  study,  symmetric  flat  composite  sandwich  panels  comprised  of  plain-weave  carbon  fabric 
preimpregnated  in  epoxy  resin  (NEWPORT  NB321/3K70P)  facesheets  and  Plascore  Nomex 
honeycomb  (PN2-3/16-3.0/4. 5/6.0)  cores  (test  section  dimensions,  8.0  in  x  8.0  in)  with  clamped 
edges  were  subjected  to  drop-weight  vertical  impact  with  a  spherical  steel  impactor  followed  by 
CAI  testing  [5].  Three  different  facesheet  configurations,  (X\  =  2  plies  [90/45]],  4  plies 
[90/45]2,  6  plies  [90/45]3),  core  densities  (X2  =  3.0,  4.5,  6.0  lb/ft3),  and  core  thicknesses  (X3  = 
3/8,  3/4,  and  9/8  in)  were  considered  in  this  examination.  The  ranges  of  material  system 
variables  are  typical  of  those  found  in  common  sandwich  panel  applications.  In  addition,  three 
different  impact  energies  (X4  =  90.0,  120.0,  150.0  in-lbs),  impactor  diameters  (X5  =  1.0,  2.0, 
3.0  in),  and  impact  velocities  (X(,  =  65.21,  96.3,  127.39  in/s)  were  considered  in  this  study.  The 
impact  parameters  correspond  to  typical  low-velocity  impacts  associated  with  relatively  blunt 
objects.  Note  that  each  sandwich  configuration  and  impact  variable,  X{  (/=  1, 2,..,  6),  assumes 
low,  midrange  or  center  point,  and  high  values  which  are  collectively  referred  to  as  the  natural 
values  of  the  independent  variables  used  in  this  study.  Table  1  summarizes  the  range  of  material 
and  impact  parameters  considered  in  this  effort.  In  addition,  the  facesheet  thickness  may  be 
characterized  in  terms  of  the  number  of  plies  (2,  4,  6)  associated  with  each  facesheet 
configuration. 

TABLE  1.  SANDWICH  CONFIGURATION  AND  IMPACT  PARAMETERS 


. - , — — - 

Natural  Values 

Material 

Variables 

Number  of  Facesheet 
Plies,  X] 

2  [90/45], 

4  [90/45]2 

6  [90/45]3 

Core  Density, 

X2  (lb/ft3) 

3.0 

4.5 

6.0 

Core  Thickness, 

Xi  (in) 

0.375 

0.750 

1.125 

Impact  Variables 

Impact  Energy, 

X4  (in-lbs) 

90.0 

120.0 

150.0 

Impactor  Diameter, 

*5  (in) 

1.0 

2.0 

3.0 

Impact  Velocity, 

X6  (in/s) 

65.21 

96.30 

127.39 

Earlier  efforts  [5  and  13]  used  through  transmission  ultrasonic  (TTU)  C-scan  measurements  and 
maximum  residual  facesheet  indentation  measurements  to  characterize  the  degree  of  impact 
damage  induced  in  the  sandwich  panels.  TTU  C-scan  measurements  provide  a  two-dimensional 
(2D)  image  of  the  projected  area  of  the  impact-damaged  region,  as  shown  in  figure  2.  The 
damage  region  may  be  characterized  either  in  terms  of  the  area  of  the  C-scan  image  or  in  terms  of 
the  diameter  of  the  damaged  region,  D,  measured  normal  to  the  direction  of  the  applied  CAI  load. 
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The  latter  represents  an  estimate  of  the  region  around  which  stress  redistribution  due  to  the 
presence  of  damage  must  occur. 


FIGURE  2.  TTU  C-SCAN  IMAGE 

Destructive  sectioning  of  impact-damaged  sandwich  composites  comprised  of  low-density 
Nomex  honeycomb  cores  (X2  =  3.0  lb/ft3)  of  variable  thickness  (X3  =  3/8,  3/4  in)  shows  that,  for 
the  given  facesheet  configurations,  the  planar  damage  region  based  upon  TTU  C-scan 
measurements  closely  corresponds  to  the  region  over  which  Nomex  honeycomb  cell  wall 
buckling/ffacture  occurs  [5].  Furthermore,  Tomblin  et  al.  [5]  demonstrated  the  CAI  residual 
strength  of  such  sandwich  composites  correlated  reasonably  well  with  the  projected  damage  area 
based  upon  TTU  C-scan  measurements  for  normal  impacts  involving  a  range  of  impactor 
diameters  and  impact  energies.  For  example,  figure  3  shows  the  normalized  CAI  residual 
strength  as  a  function  of  projected  damage  area  for  symmetric  sandwich  composites  comprised  of 
woven  fabric  carbon-epoxy  facesheets  (X  =  2  plies  [90/45]  i,  4  plies  [90/45]2,  6  plies  [90/45]3) 
and  low-density  Nomex  honeycomb  cores  (X2  =  3.0  lb/ft3)  of  variable  thickness  (X3  =  3/8,  3/4  in) 
from  [5].  The  panels  were  impacted  with  either  1.0  or  3.0  in  outside  diameter  (OD)  impactors 
over  a  range  of  impact  energies.  The  CAI  residual  strength  of  each  damaged  panel  was 
normalized  by  the  CAI  residual  strength  of  an  undamaged  (virgin)  panel  of  identical  facesheet 
configuration.  It  is  clear  from  the  figure  that  the  normalized  CAI  residual  strength  is  a  generally 
decreasing  function  of  the  projected  damage  area,  but  with  substantial  scatter.  In  addition,  the 
greatest  degradation  in  CAI  residual  strength  is  typically  associated  with  those  panels  impacted 
with  a  3.0  in  OD  impactor.  This  latter  result  is  of  particular  concern  because  sandwich  panels 
impacted  with  relatively  blunt  impactors  often  display  relatively  low  levels  of  visible  facesheet 
damage  [5  and  13].  For  the  sandwich  panels  considered  in  this  study,  Lacy  et  al.  [13]  previously 
developed  response  surfaces  for  estimating  the  planar  dimension  of  the  region  with  crushed  core 
as  well  as  the  residual  facesheet  indentation.  For  a  detailed  description  of  the  experimental  setup 
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and  procedure,  impact  testing,  nondestructive  and  destructive  damage  evaluation,  and  CAI 
testing  associated  with  this  effort,  please  see  reference  5. 


FIGURE  3.  NORMALIZED  CAI  RESIDUAL  STRENGTH  AS  A  FUNCTION  OF  TTU 

C-SCAN  AREA  FROM  REFERENCE  5 

One  primary  goal  of  this  study  is  to  develop  empirically  based  models  (i.e.,  response  surfaces) 
that  isolate  the  influence  of  key  sandwich  configuration  parameters  (e.g.,  number  of  facesheet 
plies  or  facesheet  thickness,  X\\  core  density,  Xj,  and  core  thickness,  X-$)  and  impact  parameters 
(e.g.,  impact  energy,  X4;  spherical  impactor  diameter,  X5;  and  impact  velocity,  X6)  on  the  damage 
tolerance  properties  of  sandwich  composites.  The  diameter,  D,  of  the  TTU  C-scan  area  will 
primarily  be  used  to  characterize  the  level  of  impact  damage  (cf.,  figure  2).  Lacy  et  al.  [13] 
demonstrated  that  there  may  be  no  correlation  between  the  residual  indentation  depth  (e.g., 
visible  damage)  and  the  size  of  the  internal  core  crush  region.  Combinations  of  independent 
variables,  X,  through  X6,  leading  to  a  reduced  value  of  the  diameter,  D,  of  the  projected  C-scan 
area  for  a  given  impact  event  may  minimize  the  degradation  in  the  normalized  CAI  residual 
strength  for  a  given  panel  configuration. 

The  damage-induced  residual  strength  degradation  in  sandwich  composites  due  to  foreign  object 
impact  is  highly  configuration  and  impact  parameter  dependent.  Three  distinct  sets  of  15 
experiments  were  conducted  to  isolate  the  effects  of  ( 1 )  sandwich  configuration/layup  parameters 
(Xi  through  Xi)\  (2)  the  number  of  facesheet  plies  (X,),  impact  energy  (X4),  and  impactor 
diameter  (Xs);  and  (3)  the  number  of  facesheet  plies  (X\),  impact  energy  (X4),  and  impact  velocity 
(X(>)  on  the  damage  tolerance  characteristics  of  sandwich  composite  panels.  In  each  set  of 
experiments,  three  of  the  independent  variables  were  tested  at  the  low,  center  point,  and  high 
natural  levels  defined  in  table  1 ,  while  the  remaining  independent  variables  were  held  fixed. 
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2.1  INFLUENCE  OF  SANDWICH  CONFIGURATION  PARAMETERS  ON  THE  IMPACT 
DAMAGE  TOLERANCE  OF  SANDWICH  COMPOSITES. 


Similar  to  reference  13,  a  3k  fractional  factorial  design  of  experiments  (DOE)  approach  [14]  was 
used  to  examine  the  nonlinear  interaction  effects  between  relevant  sandwich  panel  design 
parameters  and  their  influence  on  the  CAI  residual  strength,  Nyy,  associated  with  a  given  impact 
event.  To  isolate  the  coupled  effects  of  sandwich  configuration  parameters  (i.e.,  number  of 
facesheet  plies,  X\,  core  density,  Xi\  core  thickness,  X3)  on  the  damage  tolerance  characteristics 
of  sandwich  composite  panels,  the  impact  energy  (X4  =  120.0  in-lbs),  impactor  diameter 
(X5  =  3.0  in),  and  impact  velocity  (X6  =  96.30  in/s)  were  held  fixed  in  this  examination.  Based  on 
a  number  of  carefully  selected  experiments,  statistically  reliable  polynomial  expressions 
characterizing  the  impact  damage  response  of  sandwich  composites  may  be  determined  as  a 
continuous  function  of  relevant  test  parameters  [cf.,  14-19].  A  quadratic  response  surface 
generally  will  be  of  the  form 


Y  =  K  +  Z ' bt  ■  xt  +ZXV ■ xv  W 

1=1  z= 1  7=1 

where  Y  is  the  predicted  response  quantity  of  interest  (e.g.,  residual  strength);  x,  are  continuous 
normalized  independent  variables;  b0,  bh  and  by  are  least  squares  regression  coefficients;  and  k 
denotes  the  number  of  independent  variables  considered  in  the  test.  The  three  natural  values  of 
the  independent  variables,  X,  (corresponding  to  low,  center  point,  and  high  values)  are  mapped 
into  three  nondimensionalized  or  coded  levels,  x„  corresponding  to  (-1,  0,  1),  respectively.  In 
general,  the  relationship  between  the  natural  value  and  coded  level  of  a  given  independent 
variable  may  be  expressed  as 


_  2  •  Xj  (XjMAX  -YjMfN )  (2) 

(XiMAx  —XiMIN) 

where  Xnvm  and  Xmax  correspond  to  the  low  and  high  natural  values  of  the  z'th  independent 
variable,  respectively.  Equation  1  represents  a  quadratic  surface  in  the  ^-dimensional  space  of 
input  parameters;  (xb  x2,  ...  x„  ...,  xk)  =  (0,  0,  ...,  0,  ...,  0)  defines  the  center  point  of  the  test 
matrix  design  where  the  independent  variables  are  tested  at  their  midrange  values.  The  constant 
bo  represents  the  average  response  at  the  center  point  of  the  design,  while  bi  and  by  contribute  to 
the  deviation  from  this  average  value  at  points  removed  from  the  center  point  (i,  j  =1,2,...,  k). 
The  test  matrices  are  selected  so  the  generated  response  surfaces  exhibit  a  high  degree  of 
orthogonality  (i.e.,  only  bo  and  bu  are  correlated)  and  rotatability  (i.e.,  the  variance  of  the 
predicted  values  at  all  points  equidistant  from  the  center  point  is  constant).  One  key  advantage  of 
this  approach  is  that  it  can  estimate  the  impact  damage  responses  associated  with  specific 
sandwich  panel  configurations,  which  was  not  originally  tested.  In  addition,  the  coupling  and 

interactions  between  input  parameters  and  their  influence  on  the  desired  response,  Y ,  may  be 
directly  inferred  from  the  magnitudes  of  the  coefficients,  by  (i  *j). 
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The  test  matrix  developed  for  this  investigation  is  based  on  the  Box-Behnken  fractional  factorial 
design  of  experiments  technique  which  uses  the  minimum  number  of  tests  required  to  generate  a 
second-order,  statistically  reliable,  polynomial  expression  characterizing  the  response  function  of 
interest  (e.g.,  residual  strength)  [14].  Similar  approaches  have  been  used  to  characterize  the 
fracture  toughness,  residual  strength,  and  damage  tolerance  characteristics  of  composite 
laminates  [20-22].  Table  2  summarizes  the  natural  values  and  coded  levels  of  the  (k  =  3) 
independent  variables  considered  in  this  study.  Defining  the  CAI  residual  strength  as  the  desired 
response  quantity,  Nw ,  equation  1  may  be  expressed  in  expanded  form,  i.e., 

ftyy=b0+b\-X\+  b2  X2  +  W  *3  +  bi  f  Xl  +  b22  X2  +  b33  Xl 
... + bl2-  x,-  x2  +  bl2-  x,-  x3  +  b22  x2  x3 

where  x,  is  the  coded  level  of  the  ixh  independent  variable  (i  =  1,  2,  3).  Note  that  the  estimate  of 
residual  strength  in  equation  3  is  a  continuous  function  of  independent  variables.  Such  a 
representation  may  be  inappropriate  where  the  range  of  independent  variables  spans  a  bifurcation 
between  failure  modes  (e.g.,  a  transition  from  facesheet  fracture  to  facesheet  dimple 
propagation).  This  may  be  a  serious  concern  where  changes  in  facesheet  penetration  resistance 
associated  with  variations  in  the  number  of  facesheet  plies,  impactor  diameter,  or  impact  energy 
may  lead  to  multiple  CAI  failure  mechanisms. 

TABLE  2.  NATURAL  VALUES  AND  CORRESPONDING  CODED  LEVELS  OF 

MATERIAL  SYSTEM  VARIABLES 


i 

Material  Variable 

Natural  Value,  Xi 

Coded  Level,  x, 

Number  of 
Facesheet  Plies 

2  F90/451 

i 

-1 

i 

4  [90/451 

2 

0 

6  [90/451 

3 

+1 

3.0  lb/ft 

-1 

2 

Core  Density 

4.5  lb /ft* 

0 

6.0  lb/ft3 

+1 

0.375  in 

-1 

3 

Core  Thickness 

0.750  in 

0 

1.125  in 

+1 

Consistent  with  reference  13,  the  experimental  matrix  used  in  this  study  to  solve  for  the 
regression  coefficients  (b0,  bi,  and  bi/,  i,j  =  1,  2,  3)  required  that  the  high-  and  low-coded  levels 
of  any  two  independent  variables  be  paired  in  all  possible  combinations  while  fixing  the  third 
independent  variable  at  its  coded  center  point  level.  These  tests  are  performed  in  combination 
with  three  additional  experiments  in  which  the  independent  variables  are  fixed  at  their  center 
point  values  (xi  =X2  =X3  =  0),  resulting  in  a  total  of  15  experiments.  The  latter  three  center  runs 
are  important  for  assessing  the  degree  of  curvature  in  the  response  as  well  as  for  evaluating 
model  error  and  goodness  of  fit.  Table  3  shows  the  test  matrix  used  in  this  evaluation  in  terms  of 
the  coded  levels  of  the  independent  variables  as  well  as  the  measured  residual  strength,  N)y,  - 
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associated  with  each  test.  The  regression  coefficients  should  be  selected  so  the  error  between  the 
predicted  values  of  residual  strength,  Nyy,  and  the  experimentally  observed  values,  Nm  is 
minimized  over  the  range  of  panel  configurations  tested. 

Using  the  methodology  outlined  by  Lacy  et  al.  [13],  the  vector  of  desired  regression  coefficients, 

{^}=  (^0  )  ^1  ’  ^2  ’  ^3  ’  ^1 1  ’  ^22  -  ^33  ’  ^1 2  ’  3  ’  ^23  )  (^) 


may  be  determined  from 


(5) 


where  {a^,  }=  (Nyy] ,  Nvy2  ,...,NyyU  )T  is  a  vector  containing  the  experimentally  observed 

responses  (i.e.,  the  measured  residual  strengths  from  the  fifth  column  in  table  3),  and  the  matrix, 
x,  has  components 


[x]= 


t  O)  (1) 

1  x,  x2 


(*,'")• 

(*“)’ 

<0 

(1) 

(I) 

(1) 

(1) 

0) 

+  ' 

*2 

*,  ■ 

*3 

*2  ' 

*3 

(*")■ 

(*”)' 

(2) 

(2) 

(2) 

(2) 

(2) 

(2) 

*1  ' 

*1  ’ 

*3 

*3 

(  (I5)V 

(15) 

(15) 

(15) 

(15) 

(15) 

(15) 

I**  ) 

*2 

*1  ’ 

*3 

*2  ' 

*3 

(6) 


where  xfk)  represents  the  specified  coded  level  of  the  ith  independent  variable  (i  =  1,  2,  3)  during 
the  kth  experimental  run  (k  =  1, 2, . . 15)  from  table  3. 


The  vector,  b,  is  commonly  referred  to  as  the  least  squares  estimate  of  the  model  parameters 
vector.  This  estimate  is  unbiased  and  has  the  minimum  variance  among  all  other  types  of 
unbiased  estimators  [17].  Using  equations  5  and  6  in  combination  with  the  vector,  Nyy,  of 
observed  responses  from  table  3,  the  quadratic  response  surface  in  equation  3  characterizing  the 
CAI  residual  strength  as  a  continuous  function  of  coded  levels  of  the  composite  sandwich  panel 
facesheet  thickness  (xi),  core  density  fo),  and  core  thickness  (X3)  may  be  expressed  as 

Nyy  =  2149  +  937. 0-x,  +103.0-x2  +424.9-x3 

...+  5.770  •  x,2  — 190.3  •  x2  -  265.3  •  x3  (lbs/in)  (7) 

...- 55.09  x,  •  x2  +  442.1  •  x,  •  x3  +  6.380  •  x2  •  x3 


The  constant  term  (bo  =  2149  lbs/in)  represents  the  mean  CAI  residual  strength  for  those  panels 
corresponding  to  the  center  point  of  the  design.  Note  that  based  upon  the  magnitudes  of  the 
coefficients  of  the  linear  terms  in  equation  7,  the  coded  number  of  facesheet  plies,  xj,  contributes 
the  most  to  the  linear  variation  in  the  predicted  response  at  points  removed  from  the  center  point 
(xi,X2,  X3)  =  (0,  0, 0)  of  the  designed  test  matrix  (e.g.,  increasing  the  number  of  plies  from  its 
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midrange  value  will  result  in  a  linear  increase  in  the  predicted  residual  strength  from  60; 
conversely,  decreasing  the  number  of  plies  from  its  midrange  value  will  result  in  a  linear 
decrease  in  the  predicted  response).  It  makes  sense  that,  for  a  given  impact  event,  panels  with  a 
greater  number  of  facesheet  plies  will  have  a  relatively  higher  CAI  residual  strength.  Similarly, 
an  increase  in  the  core  density,  x2,  or  in  the  core  thickness,  x3,  will  contribute  to  a  linear  increase 
in  the  residual  strength.  Based  upon  the  coefficients  of  the  quadratic  terms,  variations  in  the 
coded  core  density  or  core  thickness  away  from  their  midrange  values  will  result  in  a  modest 
quadratic  decrease  in  the  predicted  CAI  residual  strength. 

Due  to  the  normalization  procedure  used  in  defining  the  coded  variables  (equation  2),  the 
influence  of  the  quadratic  terms  on  the  predicted  response  will  be  maximized  as  the  coded 
variables  approach  their  extreme  values  (*,  =  ±1).  The  coefficients  of  the  interaction  terms 
(x,  •  x/,  i  *  j)  provide  an  indication  of  the  complex  coupling  between  material  parameters  and 
their  influence  on  the  sandwich  panel  damage  tolerance  characteristics.  In  this  case,  a 
simultaneous  increase  in  the  number  of  facesheet  plies  and  core  thickness  will  have  a  beneficial 
effect  on  the  predicted  response.  Similar  arguments  can  be  made  when  interpreting  the  influence 
the  remaining  terms  in  equation  7  has  on  the  predicted  response.  It  is  clear  from  the  functional 
form  of  the  response  surface,  (equation  7),  that  for  the  given  range  of  material  and  impact 
parameters,  the  CAI  residual  strength  is  somewhat  sandwich  configuration  dependent. 

Note  that  the  Box-Behnken  test  matrix  defined  in  table  3  represents  a  spherical  experimental 
design,  i.e.,  the  first  12  tests  are  conducted  at  points  (xi,  x2,  x3)  that  are  equidistant  from  the 
center  point  of  the  design.  If  one  considers  a  cube  in  the  space  of  the  coded  variables  defined  by 
the  planes  (x}  =±1,  x2  =  ±1,  x3  =  ±1),  the  noncenter  point  tests  are  conducted  at  edge  points  (i.e., 
points  located  at  the  centers  of  the  edges  of  the  cube)  as  shown  schematically  in  figure  4.  Hence, 

the  edge  points  all  lie  a  coded  radial  distance,  r  =  Jxf  +  x22  +  x2  =S  =  1.414,  from  the  center 

point  of  the  design.  When  using  response  surfaces,  such  as  equation  7,  to  estimate  the  impact 
damage  response  associated  with  panel  configurations  not  considered  in  the  original  test  matrix, 
special  care  must  be  used  to  ensure  that  the  interpolation  occurs  within  the  spherical  domain 

defined  by  the  coded  radius,  r  =  yjl .  Extrapolation  to  panel  configurations  that  lie  outside  of 
this  spherical  region  may  lead  to  serious  errors  in  the  predicted  response.  If  response  surface 
estimates  are  desired  over  a  cuboidal  domain,  then  various  other  experimental  designs  may  be 
used  (e.g.,  Central  Composites  Design;  cf.  [17]). 

While  the  interpretation  of  the  regression  coefficients  is  most  easily  performed  if  a  response 
surface  is  characterized  in  terms  of  the  coded  levels  of  the  independent  variables,  equation  2  may 
be  used  to  rewrite  the  response  surface  in  equation  7  in  terms  of  the  natural  values  of  the 
independent  variables,  i.e., 

Nyy  =-2159  +  97.56  X,  +894.7 -X2  + 1554- JST3 

-  + 1  -44 1  •  X2  -  84.56  •  X2  - 1 887  •  X2  (lbs/m)  (8) 

...-18.36-JT,  X2  +589.4  X,  X3  + 11.35  -X2  -X3 


8 


The  latter  representation  allows  the  response  surface  to  be  plotted  in  the  space  of  the  natural 
variables.  The  response  surface  (equation  7  or  8)  may  be  used  to  identify  those  panel 
configurations  that  retain  a  higher  fraction  of  their  virgin  panel  strength  for  a  given  impact  event. 


Table  3  summarizes  the  measured  and  predicted  CAI  residual  strength  for  the  sandwich 
composite  panels  considered  in  this  study.  The  difference  between  the  experimentally  measured 
CAI  strengths  and  the  estimated  values  using  equation  7  varied  between  0.4%  and  9.1%  for  the 
15  panels  tested,  with  a  mean  difference  of  3.4%.  This  suggests  that  the  response  surface 
(equation  7)  may  provide  a  reasonable  characterization  of  the  influence  of  the  number  facesheet 
plies,  core  density,  and  core  thickness  on  CAI  residual  strength.  It  is  useful  to  recognize  that 
those  specimens  with  two-ply  facesheet  configurations  (xi  =  -1)  showed  the  biggest  disparity 
between  the  measured  and  estimated  CAI  residual  strengths.  The  regression  results  may  also  be 
viewed  schematically  in  the  space  of  coded  independent  variables,  as  shown  in  figure  4.  Note 
that  the  experimental  center  point  observations  (cfi,  table  3)  are  tightly  banded  about  their  mean 
value,  bo  =  2149  lbs/in.  This  is  highly  desirable  in  order  to  avoid  introducing  a  large  degree  of 
uncertainty  into  the  regression  model.  This  can  be  a  particular  concern  where  the  experimental 
observations  may  not  be  particularly  repeatable.  Standard  analysis  of  variance,  error  estimation, 
and  goodness  of  fit  checks  [17]  suggest  that  the  response  surface  in  equation  7  provides  a 
statistically  reliable  estimate  of  CAI  residual  strength. 

TABLE  3.  CODED  SANDWICH  CONFIGURATION  VARIABLES  TEST  MATRIX  AND 
MEASURED  AND  PREDICTED  RESIDUAL  STRENGTH 


Test, 

k 

X\ 

X2 

*3 

Measured 
Residual 
Strength, 
Nyy  (lbs/in) 

Predicted 

Residual 

Strength, 

Nyy  (lbs/in) 

Nyy 

1%) 

1 

+1 

+i 

0 

2996 

2949 

1.6 

2 

+1 

-l 

0 

2952 

2853 

3.3 

3 

-1 

+i 

0 

1087 

1185 

9.1 

4 

-1 

-l 

0 

822.2 

868.9 

5.7 

5 

+1 

0 

+1 

3637 

3693 

1.5 

6 

+1 

0 

-1 

1870 

1959 

°o 

7 

-1 

0 

+1 

1024 

934.7 

8.7 

8 

-1 

0 

-1 

1025 

969.2 

5.5 

9 

0 

+1 

+1 

2236 

2227 

0.4 

10 

0 

+1 

-1 

1407 

1365 

3.0 

11 

0 

-1 

+1 

1966 

2008 

2.2 

12 

0 

-1 

-1 

1162 

1172 

0.8 

13 

0 

0 

0 

2107 

2149 

2.0 

14 

0 

0 

0 

2182 

2149 

1.5 

15 

0 

0 

0 

2157 

2149 

0.4 

o 


1407 

1365 


1024 

934.7 


1870 

1959 


*2  -  Coded  core 
*  density 


=  Coded  number  of 
facesheet  plies 


X3  =  Coded  core 
thickness 


Nyy=  Measured  residual  strength  (Ibs/ln) 
A^yy=  Estimated  residual  strength  (Ibs/ln) 


FIGURE  4.  RESPONSE  SURFACE  ESTIMATES  OF  RESIDUAL  STRENGTH  AS  A 
FUNCTION  OF  SANDWICH  CONFIGURATION  PARAMETERS 

Figure  5  summarizes  the  influence  the  number  of  facesheet  plies  (Aj,  x,),  core  density  (X2,x2), 
and  core  thickness  (X3,  x3)  have  on  the  estimate  of  CAI  residual  strength  (equation  8)  for  a  fixed 
impact  energy,  impactor  diameter  and  velocity.  In  the  response  surface  plots,  the  bounds  of  the 

sphere  of  coded  radius,  r  =  yfl,  is  denoted  by  an  inscribed  dashed  circle.  Strictly  speaking, 
regression  data  falling  outside  the  bounds  of  this  sphere  would  correspond  to  an  extrapolation  of 
the  response  surface  results.  Figures  5(a),  5(b),  and  5(c)  show  the  CAI  residual  strength 
(equation  8)  as  a  function  of  core  density  and  core  thickness  for  the  two-ply,  four-ply,  and  six-ply 
facesheet  configurations  (or  x\  =  -1,  0,  1),  respectively.  A  comparison  of  the  three  figures  reveals 
that  the  predicted  strength  increases  substantially  as  the  number  of  facesheet  plies  is  increased.  It 
is  clear  from  figure  5(a),  that  for  sandwich  panels  with  the  minimum  number  of  facesheet  plies 
{X\  =  2  plies;  xj  =  -1),  the  CAI  strength  generally  increases  with  increasing  core  density,  X2.  The 
CAI  strength  is  a  relative  maximum  in  the  vicinity  of  the  maximum  core  density  (X2  =  6.0  lb/ft3; 
Xl  7  P  anc^  Grange  core  thickness  (X$  =3/4  in;  x3  =  0).  The  minimum  strength  occurs  in  the 
vicinity  of  X2  —  3.0  lb/ft  (x2  =  -1)  and  X2  =  3/4  in  (x3  =  0).  The  relatively  reduced  penetration 
resistance  associated  with  the  two-ply  facesheet  configuration  likely  results  in  more  localized 
damage  development  for  a  given  impact  event,  particularly  as  the  core  density  is  decreased. 

As  the  number  of  facesheet  plies  is  increased  to  its  midrange  value  (X\  =  4  plies;  x,  =  0),  the 
magnitude  of  the  peak  CAI  strength  and  its  location  in  the  space  of  independent  variables 
changes  (figure  5(b)).  The  estimated  CAI  strength  tends  to  increase  rapidly  with  increasing  core 
thickness  and  is  less  sensitive  to  the  density  of  the  core.  For  the  relatively  blunt  object  impacts 
considered  here,  enhanced  penetration  resistance  associated  with  an  increase  in  the  number  of 
facesheet  plies  often  results  in  CAI  failure  processes  characterized  by  facesheet  dimple 
propagation  across  the  width  of  the  test  specimen  [5].  The  applied  load  level  (e.g.,  strength),  at 
which  this  process  initiates,  apparently  is  influenced  by  core  thickness  (e.g.,  increasing  the  core 
thickness  may  provide  enhanced  bending  rigidity  that  serves  to  delay  the  onset  of  facesheet 
dimple  propagation).  Here  the  peak  predicted  CAI  strength  occurs  in  the  vicinity  of 
X2  =  6.0  lb/ft  (x2=  1)  and  X2  =  9/8  in  (x3  =  1);  the  strength  is  a  relative  minimum  for 
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X2  =  3.0  lb/ft3  (x2=  -1)  and  X3  =  3/8  in  (jc3  =  -1),  corresponding  to  the  extreme  ranges  of  the 
tested  values. 
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FIGURE  5.  PREDICTED  RESIDUAL  STRENGTH  (lbs/in):  (a)  Xx  =  2  plies  [90/45] i, 
(b)  Xx  =  4  plies  [90/45]2,  AND  (c)  Xi  =  6  plies  [90/45]3 
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Figure  5(c)  shows  that  for  sandwich  panels  with  six-ply  facesheet  configurations  (Xt  =  6  plies; 

=  1),  the  estimated  CAI  strength  increases  significantly  with  increasing  core  thickness  and  is 
fairly  insensitive  to  core  density.  There  is  likely  a  competition  between  enhanced  penetration 
resistance  and  improved  bending  stiffness  associated  with  changes  in  independent  variables  that 
govern  the  damage  development  [  1 3]  and  residual  strength. 

The  preceding  results  are  in  agreement  with  the  experimental  observations.  Figures  6  through  8 
show  the  recorded  CAI  residual  strengths  for  the  test  specimens  corresponding  to  the  midrange 
core  thickness  configurations  (*3  =  0;  X$  =  3/4  in),  midrange  core  density  configurations  (*2  =  0; 
X2  =  4.5  lb/ft3),  and  midrange  facesheet  configurations  (*i  =  0;  X\  =  4  plies),  respectively. 
Included  with  each  observation  is  its  associated  TTU  C-scan  image  and  measured  diameter,  D, 
from  reference  13.  These  figures  suggest  that  the  experimentally  observed  CAI  strengths  for  a 
given  facesheet  configuration  are  consistent  with  the  regressions  in  figure  5.  Furthermore,  the 
observed  CAI  strengths  for  the  four-ply  and  six-ply  facesheet  configurations  generally  were 
inversely  proportional  to  the  diameter  of  the  TTU  C-scan  images  reported  in  reference  13. 

For  the  test  matrix  summarized  in  table  3,  figure  9  contains  a  comparison  of  the  response  surface 
predictions  of  the  TTU  C-scan  diameter  from  reference  13  to  the  response  surface  predictions  of 
CAI  residual  strength  developed  in  this  examination.  For  illustration  purposes,  the  CAI  residual 
strength  of  each  damaged  panel  was  normalized  by  the  CAI  residual  strength  of  an  undamaged 
(virgin)  panel  of  identical  facesheet  configuration  prior  to  determination  of  the  regression 
coefficients  (equation  5).  A  comparison  of  figures  9(a),  9(b),  and  9(c)  to  figures  9(d),  9(e),  and 
9(f),  respectively,  suggests  that  with  the  exception  of  the  two-ply  facesheet  configuration  (figures 
9(a)  and  9(d)),  the  predicted  normalized  CAI  residual  strength  is  generally  inversely  proportional 
to  the  size  of  the  planar  damage  area,  D,  for  the  given  impact.  This  is  consistent  with  the 
observations  of  Tomblin  et  al.  [5]  (cf.,  figure  2).  Lacy  et  al.  [13]  noted  that  the  reduced 
penetration  resistance  associated  with  the  two-ply  facesheet  configurations  likely  resulted  in  a 
decrease  in  the  planar  size  of  the  internal  damage  for  a  given  impact  event.  There  may  be, 
however,  appreciably  more  localized  facesheet  damage  (ply  fracture,  fiber  breaks,  matrix  cracks,' 
etc.)  that  reduces  the  facesheet  load  transfer  across  the  impact  site.  The  loss  of  load  carrying 
capability  of  the  facesheet  coupled  with  the  reduced  bending  stiffness  associated  with  the  two-ply 
facesheet  configuration  results  in  a  bifurcation  between  compressive  failure  modes 
(i.e.,  facesheet  fracture  occurs  rather  than  progressive  facesheet  dimple  propagation).  Hence, 
there  may  be  no  clear  correlation  between  the  measured  TTU  C-scan  diameter  and  residual' 
strength  for  impacts  involving  partial  facesheet  penetration.  Clearly,  using  continuous  function 
response  surfaces  to  predict  residual  strength  should  be  used  with  extreme  caution  over  ranges  of 
independent  variables  where  a  potential  bifurcation  between  failure  modes  may  result  in  an 
apparent  step  discontinuity  in  the  actual  residual  strength,  particularly  for  those  impacts 
involving  smaller  diameter  impactors,  increasing  impact  energy,  and/or  thinner  facesheets.  This 
may  explain  why  the  disparity  between  the  measured  and  predicted  CAI  residual  strengths  was 
the  greatest  for  sandwich  composites  with  two-ply  facesheet  configurations  (cf.,  table  3).  Similar 
to  the  results  shown  in  figure  5,  the  estimated  residual  strength  for  the  four-ply  and  six-ply 
facesheet  configurations  may  increase  significantly  with  increasing  core  thickness. 
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FIGURE  6.  MEASURED  DAMAGE  DIAMETER  [13]  AND  RESIDUAL  STRENGTH 
FOR  CORE  THICKNESS,  A3  =  3/4  in  (x3  =  0) 
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X2  =  Core  density  (lb/ft3) 


FIGURE  9.  PREDICTED  DAMAGE  DIAMETER  (in)  FOR Xx:  (a)  TWO,  (b)  FOUR,  AND 
(c)  SIX  PLIES  [13];  PREDICTED  NORMALIZED  RESIDUAL  STRENGTH 
FOR^:  (d)  TWO,  (e)  FOUR,  AND  (f)  SIX  PLIES 
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A  comparison  of  figures  9(d),  9(e),  and  9(f)  suggests  that  the  given  impact  can  result  in  residual 
strength  values  that  range  roughly  between  50%-90%  of  the  virgin  panel  strength  for  a  given 
facesheet  configuration.  Response  surfaces  similar  to  equation  8  may  allow  manufacturers  to 
estimate  the  peak  degradation  in  CAI  residual  strength  for  a  given  class  of  impact  (i.e.,  establish 
conservative  knockdown  factors  to  the  virgin  panel  strength  that  arguably  account  the  influence 
of  expected  impacts).  Such  information  may  prove  valuable  in  the  initial  design  and  sizing  of 
key  sandwich  composite  structural  components  as  well  as  in  the  specification  of  design 
allowables  for  damage  tolerant  parts.  Special  care  must  be  taken  to  ensure  that  the  range  of 
impact  variables  (impactor  diameter,  impact  energy,  impact  velocity,  etc.)  used  in  the 
development  represent  a  realistic  threat  to  a  given  structural  element. 

Table  4  contains  a  limited  set  of  additional  independent  experimental  results  from  Tomblin  et  al. 
[5],  where  the  combination  of  sandwich  configuration  variables  (X\,  X2,  A3)  =  (4  plies,  3.0  lb/ft3, 
3/4  in)  or  (x\,  x2,  *3)  =  (0,  -1,  0)  lie  within  a  coded  radius,  r  =  yfl  =1.414 ,  of  the  center  point  of 

the  experimental  design.  Response  surface  predictions  of  the  planar  diameter,  D ,  of  the  internal 
damage  from  reference  13  are  compared  to  the  measured  diameter,  D,  and  the  CAI  residual 
strength  (equation  7)  are  compared  to  measured  residual  strengths  in  table  4.  Note  that  response 
surface  estimates  for  predicted  CAI  residual  strength  using  equation  7  more  closely  matched  the 
limited  experimental  results  than  response  surface  estimates  of  the  TTU  C-scan  diameter  from 
reference  13.  This  is  consistent  with  scatter  observed  in  figure  3.  The  mean  value  of  the 
experimentally  CAI  residual  strengths  {Nyymean  =  1792  lbs/m)  for  this  sandwich  panel 
configuration  differed  from  the  value  predicted  in  equation  7  by  3.5%.  The  difference  between 
individual  measurements  and  the  predicted  value  varied  between  3.7%  and  8.7%,  with  a  mean 
difference  of  5.9%.  This  suggests  that  response  surfaces  similar  to  equation  7  may  be  useful 
tools  in  identifying  combinations  of  sandwich  composite  parameters  leading  to  improved  damage 
tolerance  properties  for  a  given  class  of  impact. 

TABLE  4.  INTERPOLATION  OF  REGRESSION  RESULTS  IN  THE  SPACE  OF  CODED 

MATERIAL  VARIABLES 


Xi 

*2 

*3 

Measured 
Damage 
Diameter, 
D  (in) 
[13] 

Predicted 

Damage 

Diameter, 

D  (in) 
[131 

D-D 

D 

(%) 

Measured 

Residual 

Strength, 

Nyy 

(lbs/in) 

Predicted 

Residual 

Strength, 

K 

(lbs/in) 

Nyy-Nyy 

Nyy 

(%) 

0 

-1 

0 

1 

2.046 

2.239 

9.4 

1736 

1855 

6.9 

0 

-1 

0 

1 

2.131 

2.239 

5.1 

1935 

1855 

4.1 

0 

-1 

0 

1 

2.056 

2.239 

8.9 

1790 

1855 

3.7 

0 

-1 

0 

1 

2.851 

2.239 

21.5 

1708 

1855 

8.7 
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22  INFLUENCE  OF  FACESHEET  THICKNESS.  IMPACT  ENERGY.  AND  IMPACTOR 
DIAMETER  on  the  impact  damage  tolerance  OF  SANDWICH  COMPOSITES 

In  the  preceding  section,  the  isolated  influence  of  sandwich  configuration  parameters  on  the 
damage  tolerance  properties  of  sandwich  composites  was  investigated  for  a  fixed  set  of  impact 
parameters.  In  this  study,  the  isolated  effects  of  the  number  of  facesheet  plies  (A'i),  impact 
energy  (X4),  and  impactor  diameter  (V5)  on  the  damage  tolerance  characteristics  of  sandwich 
composite  panels  were  investigated  using  the  Box-Behnken  experimental  design.  Lacy  et  al.  [13] 
conducted  a  similar  study  investigating  the  influence  of  these  parameters  on  the  damage 
resistance  characteristics  of  sandwich  composites  for  the  same  set  of  test  specimens  considered 
here.  Table  5  summarizes  the  low,  midrange,  and  high  levels  of  the  natural  (Xu  X4,  X5)  and 
coded  (x,,  x4,x5)  independent  variables  considered  in  this  examination.  Here,  the  core  density 

(X2  =  3.0  lb/ft ),  core  thickness  (X3  =  3/4  in),  and  impact  velocity  (X6  =  96.3  in/s)  were  held 
fixed. 

TABLE  5.  NATURAL  VALUES  AND  CORRESPONDING  CODED  LEVELS  OF  THE 
FACESHEET  PLIES  AND  IMPACT  VARIABLES 


i 

Independent  Variable 

Natural  Value,  X,- 

Coded  Level,  x,- 

1 

Number  of  Facesheet 
Plies 

2  [90/451, 

-1 

4  r90/4512 

0 

6  [90/451, 

+1 

Impact  Energy 

90.0  in-lbs 

-1  1 

4 

120.0  in-lbs 

0 

150.0  in-lbs 

+1 

Impactor  Diameter 

1.0  in 

-i  J 

5 

2.0  in 

0 

3.0  in 

+1 

Table  6  summarizes  the  combinations  of  coded  independent  variables  as  well  as  the 
experimentally  measured  residual  strengths  for  the  15  experiments  used  in  the  regression 
analysis.  Following  a  procedure  similar  to  that  used  to  develop  the  response  surfaces  (equation  7 
and  8),  statistically  reliable,  second-order  response  surfaces  were  generated  that  characterize  the 
CAI  residual  strength  as  a  continuous  function  of  the  number  of  facesheet  plies,  impact  energy, 
and  impactor  diameter.  It  is  especially  important  to  consider  the  possibility  of  bifurcations 
between  failure  modes  when  interpreting  the  regression  results.  An  estimate  of  the  CAI  residual 
strength  from  the  regression  analysis  may  be  expressed  either  in  terms  of  coded  or  natural  values 
of  the  independent  variables,  i.e., 

Nyy  =  1 722  + 733.0 -x,  -  225.2  -x4  -35.57  x5 

...  - 1 34.0  •  x2  +  29.57  ■  x2  - 1 77. 1  •  x2  (lbs/in)  (9a) 

...  — 221.0-x,  -x4  — 161.6  x,  x5  -72.47 -x4  -x5 
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and 


N„  =-2537  +  1238-  Xx  +  4. 176-  X4+ 1286-  X5 

...-33.50-X2  +  0.033-X2  -177.1- A2  (lbs/m)  (9b) 

...-3.684  ■XI-X4-  80.78  •  Xt  ■  X,  -  2.4 1 6  •  X4  ■  X5 

TABLE  6.  COMPARISONS  BETWEEN  PREDICTED  AND  MEASURED 

RESIDUAL  STRENGTH 


Test, 

k 

*1 

X4 

x5 

Measured 
Residual 
Strength, 
Nn,  (Ibs/in) 

Predicted 

Residual 

Strength, 

N„  (lbs/in) 

N  -N 

yy  yy 

<%)  _ 

16 

+1 

+1 

0 

1913 

1905 

0.4 

17 

+1 

-1 

0 

2942 

2797 

4.9 

18 

-1 

+1 

0 

735.5 

880.7 

19.7 

19 

-1 

-1 

0 

880.5 

889 

1.0 

20 

+1 

0 

+1 

1811 

1947 

7.5 

21 

+1 

0 

-1 

2323 

2341 

0.8 

22 

-1 

0 

+1 

822.1 

804.2 

2.2 

23 

-1 

0 

-1 

688.0 

552.2 

19.7 

24 

0 

+1 

+1 

1369 

1242 

9.3 

25 

0 

+1 

-1 

1467 

1458 

0.6 

26 

0 

-1 

+1 

1827 

1837 

0.5  j 

27 

0 

-1  i 

-1 

1636 

1763 

7.8 

28 

0 

0 

0 

1666 

1722 

3.4 

29 

0 

0 

0 

1712 

1722 

0.6 

30 

0 

0 

0 

1789 

1722 

3.7 

The  constant  term  (bo  =  1 722  lbs/in)  in  equation  9a  represents  the  mean  CAI  residual  strength  for 
those  panels  corresponding  to  the  center  point  of  the  design.  Based  upon  the  magnitudes  of  the 
coefficients  of  the  linear  terms  in  equation  9a,  the  coded  number  of  facesheet  plies,  x\,  and 
impact  energy,  X4,  contribute  the  most  to  the  linear  variation  in  the  predicted  response  at  points 
removed  from  the  center  point  (x\,  X4,  x5)  =  (0,  0,  0)  of  the  designed  test  matrix.  As  might  be 
expected,  increasing  the  number  of  facesheet  plies  and/or  decreasing  the  impact  energy  and 
impactor  diameter  will  result  in  a  linear  increase  in  the  estimated  CAI  strength  from  bo .  Lacy  et 
al.  [13]  noted  that  increasing  either  the  impact  energy  or  the  impactor  diameter  will  generally 
produce  internal  damage  that  is  spread  over  a  greater  area.  Similar  arguments  can  be  made  when 
interpreting  the  influence  of  the  remaining  quadratic  and  coupling  terms  in  equation  9a  on  the 
predicted  response. 
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2.2. 1  Discussion  of  Facesheet  Plies  and  Impact  Parameter  Regression  Model  Results. 


Table  6  summarizes  the  measured  and  predicted  CAI  residual  strength  for  each  of  the  sandwich 
composite  panels  considered  in  this  effort.  The  difference  between  the  experimentally  measured 
residual  strengths  and  the  estimated  values  using  equation  9  varied  between  0.4%  and  19.7%  for 
the  15  panels  tested,  with  a  mean  difference  of  5.5%.  This  suggests  that  the  response  surface 
(equation  9)  may  provide  a  reasonable  characterization  of  the  influence  of  the  number  facesheet 
plies,  impact  energy,  and  impactor  diameter  on  the  CAI  residual  strength  if  the  number  of 
facesheet  plies  is  greater  than  two.  The  regression  results  may  also  be  viewed  schematically  in 
the  space  of  coded  independent  variables,  as  shown  in  figure  1 0.  Similar  to  the  previous  results 
(cfi,  table  3),  the  differences  between  the  observed  and  predicted  residual  strengths  were  greatest 
for  those  sandwich  panels  with  two-ply  facesheets  (-20%  disparity).  Recall  that  sandwich  panels 
comprised  of  thinner  facesheets  are  somewhat  more  susceptible  to  failure  in  compression  due  to 
facesheet  fracture,  since  the  relative  lack  of  bending  rigidity  in  the  individual  facesheets  may 
impede  the  onset  of  progressive  dimple  propagation.  The  likelihood  of  facesheet  fracture  maybe 
exacerbated  by  any  stress  concentration  due  to  varied  facesheet  damage.  As  the  number  of 
facesheet  plies  is  reduced,  and  the  possibility  of  facesheet  penetration  increases,  there  is  a 
transition  in  failure  mode  from  predominately  progressive  dimple  propagation  to  facesheet 
fracture.  This  may  explain  why  the  regression  results  do  not  more  closely  correlate  with  the 
experimental  observations  for  panels  with  two-ply  facesheets. 
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FIGURE  10.  RESPONSE  SURFACE  ESTIMATES  OF  RESIDUAL  STRENGTH  AS  A 
FUNCTION  OF  FACESHEET  PLIES  AND  IMPACT  PARAMETERS 

Figure  11  summarizes  the  influence  of  the  number  of  facesheet  plies  (Xu  *,),  impact  energy 
(X4,  x4),  and  impactor  diameter  (X5,  x5)  on  the  estimate  of  CAI  residual  strength  (equation  9). 
Figures  11(a),  11(b),  and  11(c)  show  the  CAI  residual  strength  (equation  9)  as  a  function  of 
impact  energy  and  impactor  diameter  for  the  two-ply,  four-ply,  and  six-ply  facesheet 
configurations  (or  xj  =-1,0,  1),  respectively.  Similar  to  the  results  from  the  previous  section,  a 
comparison  of  the  three  figures  reveals  that  the  predicted  strength  increases  substantially  as  the 
number  of  facesheet  plies  is  increased.  For  the  two-ply  facesheet  configuration  (figure  1 1(a)),  the 
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FIGURE  1 1 .  PREDICTED  RESIDUAL  STRENGTH  (lbs/m):  (a)  X,=2  plies  [90/45] 
(b)  X\  =  4  plies  [90/45]2,  AND  (c)  X,  =  6  plies  [90/45]3 


Increasing  number  of  facesheet  plies, 


estimated  CAI  strength  generally  increases  with  increasing  impactor  diameter  and  is  somewhat 
insensitive  to  the  impact  energy  level.  One  possible  explanation  for  the  increase  in  the  estimated 
strength  is  that  as  the  diameter  of  the  impactor  is  increased  and  the  potential  for  facesheet 
penetration  is  reduced,  more  load  may  be  transferred  across  the  impact  site.  The  compressive 
failure  mechanisms  and  residual  strength  are  likely  influenced  by  whether  or  not  partial  facesheet 
penetration  occurs.  This  is  of  particular  concern  for  those  impacts  involving  smaller  diameter 
impactors  and/or  increasing  impact  energy.  For  the  four-ply  and  six -ply  facesheet  configurations 
(figures  11(b)  and  11(c)),  which  have  somewhat  higher  penetration  resistance  and  bending 
rigidity,  the  estimated  strength  is  a  generally  decreasing  function  of  impact  energy.  For  the  four- 
ply  facesheet  configuration,  there  is  additional  damage  for  small  diameter  impactors  that  reduces 
the  residual  strength.  At  midrange-to-high  impact  energy  levels,  the  predicted  response  is  a 
decreasing  function  of  impactor  diameter,  particularly  for  the  six-ply  facesheet  configuration; 
increasing  the  impactor  diameter  from  2.0  in  to  3.0  in  OD  can  result  in  a  significant  reduction  in 
the  estimated  CAI  strength  (figure  11(c)).  For  specimens  with  thicker  facesheets,  widespread 
core  damage  resulting  from  impacts  involving  larger  diameter  impactors  may  promote  CAI 
failure  due  to  progressive  facesheet  dimple  propagation.  Small  diameter  impactors  seem  to 
introduce  other  damage  mechanisms.  This  underscores  the  need  to  consider  a  variety  of  impactor 
diameters  when  establishing  a  damage  tolerance  plan  for  sandwich  composite  aircraft  structures. 

The  preceding  results  are  in  agreement  with  the  experimental  observations.  Figures  12  through 
14  show  the  recorded  CAI  residual  strengths  for  the  test  specimens  corresponding  to  the 
midrange  impactor  diameter  (x5  =  0;  X5  =  2  in),  midrange  impact  energy  (x4  =  0;X4=  120  in-lbs), 
and  midrange  facesheet  configurations  (xi  =  0;  X\  =  4  plies),  respectively.  Included  with  each 
observation  is  its  associated  TTU  C-scan  image  and  measured  diameter,  D,  from  reference  13. 
These  figures  suggest  that  the  experimentally  observed  CAI  strengths  for  a  given  facesheet 
configuration  are  consistent  with  the  regressions  of  figure  11.  Again,  the  observed  CAI  strengths 
for  the  four-ply  and  six-ply  facesheet  configurations  generally  were  inversely  proportional  to  the 
diameter  of  the  TTU  C-scan  images  reported  in  reference  13. 

For  the  test  matrix  summarized  in  table  6,  figure  15  contains  a  comparison  of  the  response 
surface  predictions  of  the  TTU  C-scan  diameter  from  reference  13  to  the  response  surface 
predictions  of  the  normalized  CAI  residual  strength  developed  in  this  examination.  Comparing 
figures  15(a),  15(b),  and  15(c)  to  figures  15(d),  15(e),  and  15(f),  respectively,  shows  that,  with 
the  exception  of  the  two-ply  facesheet  configuration  (figures  15(a)  and  15(d)  and  four-ply 
facesheet  configurations  under  low-impact  energy  and  small  impactor  diameters  (figures  15(b) 
and  15(e)),  the  maximum  degradation  in  residual  strength  and  the  maximum  TTU  C-scan 
diameter  both  correspond  to  higher-energy  impacts  involving  3.0  in  OD  impactors.  This  is 
consistent  with  the  observations  of  Tomblin  et  al.  [5]  (cf.,  figure  2).  Again,  for  impacts 
involving  partial  facesheet  penetration,  there  may  be  no  clear  correlation  between  the  measured 
TTU  C-scan  diameter  and  CAI  residual  strength.  A  comparison  of  figures  15(d),  15(e),  and  15(f) 
suggests  that  the  CAI  residual  strength  can  range  roughly  between  45%-70%  of  the  virgin  panel 
strength  for  a  given  facesheet  configuration.  Similar  to  the  results  shown  in  figure  11,  the 
estimated  normalized  residual  strength  for  the  four-ply  and  six -ply  facesheet  configurations  may 
decrease  significantly  with  increasing  impact  energy  and  impactor  diameter. 
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FIGURE  12.  MEASURED  DAMAGE  DIAMETER  [13]  AND  RESIDUAL  STRENGTH  FOR 

IMPACTOR  DIAMETER,  X5  =  2.00  in  (x5  =  0) 
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FIGURE  13.  MEASURED  DAMAGE  DIAMETER  [13]  AND  RESIDUAL  STRENGTH  FOR 

IMPACT  ENERGY,  X4  =  1 20  in-lbs  (x4  =  0) 
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FIGURE  14.  MEASURED  DAMAGE  DIAMETER  [13]  AND  RESIDUAL  STRENGTH  FOR 
FACESHEET  CONFIGURATION,  X,  =  4  plies  [90/45]2  (x,  =  0) 
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FIGURE  15.  PREDICTED  DAMAGE  DIAMETER  (in)  FORA",:  (a)  TWO,  (b)  FOUR,  AND 
(c)  SIX  PLIES  [13];  PREDICTED  NORMALIZED  RESIDUAL  STRENGTH 
FOR  X i :  (d)  TWO,  (e)  FOUR,  AND  (f)  SIX  PLIES 
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2.2.2  Correlation  With  Independent  Experimental  Results. 


Table  7  contains  a  set  of  additional  independent  experimental  results  from  Tomblin  et  al.  [5], 
where  the  ranges  of  coded  material  system  and  impact  variables  are  in  the  vicinity  of  those 

A 

considered  in  this  study.  Response  surface  predictions  of  the  planar  diameter,  D ,  of  the  internal 
damage  from  reference  13  are  compared  to  measured  damage  diameters,  D,  from  reference  13. 
Also  shown  in  table  7  are  the  measured  and  the  predicted  (using  equation  9)  CAI  residual 
strengths.  Note  that  the  predicted  errors  for  CAI  residual  strength  using  equation  9  are 
comparable  to  those  for  the  damage  diameter  from  reference  13.  The  combinations  of  the  coded 
number  of  facesheet  plies  (*i),  impact  energy  (x4),  and  impactor  diameter  (xs),  contained  in  the 
first  eight  rows  of  the  table  lie  within  a  coded  radius,  r  =  V 2  =  1.414 ,  of  the  center  point  of  the 
experimental  design;  interpolation  of  response  surface  results  is  appropriate  for  these  data.  The 
remaining  combinations  of  independent  variables  correspond  to  r  >  a/2  .  Using  the  response 
surface  (equation  9)  for  the  latter  case  corresponds  to  an  extrapolation  beyond  the  range  of 
independent  variables  used  in  the  regression;  such  estimates  are  denoted  in  the  table  by  an 
asterisk  (*).  The  difference  between  the  experimentally  measured  CAI  residual  strengths  and  the 
interpolated  values  varied  between  1.7%  and  22.6%,  with  a  mean  difference  of  14.0%.  Hence, 
the  response  surface  estimates  using  equation  9  correlate  moderately  well  with  the  experimental 
data  within  the  space  of  the  Box-Behnken  design.  The  first  three  extrapolated  values  of  the 
strength  contained  in  the  table  correspond  to  combinations  of  independent  variables  that  fall  just 
outside  of  the  coded  radius,  r  =  a/2;  the  response  surface  estimates  (equation  9)  correlate 
reasonably  well  with  the  experimental  data  for  these  cases.  The  remaining  combinations  of 
independent  variables,  however,  correspond  to  extrapolation  of  regression  results  toward  various 

comers  of  the  cube  shown  in  figure  10  (i.e.,  r  =  y]xf  +xl+x*  ~a/3  =1.732). 

Table  7  shows  the  magnitudes  of  the  differences  between  measured  and  extrapolated  values  of 
the  CAI  residual  strength  that  can  be  significantly  larger  than  for  the  case  involving  interpolation 
of  results.  The  maximum  difference  between  the  measured  and  extrapolated  results  (46.6%) 
occurred  for  (xu  x4,  x5)  =  (-1,  -0.95,  -1)  or  (Xh  X4,  X5)  =  (2  plies,  91.5  in-lbs,  1.0  in).  Figure 
1 1(a)  shows  that  a  relatively  large  gradient  exists  in  the  extrapolated  response  surface  prediction 
along  the  radial  line  defined  by  the  points  (Xi,  X4,  X5)  =  (2  plies,  120.0  in-lbs,  2.0  in)  and  (Xi, 
X4,  X5)=  (2  plies,  90.0  in-lbs,  1.0  in).  Clearly,  extreme  caution  must  be  used  to  ensure  that 
response  surfaces  estimates  are  only  evaluated  in  the  spherical  domain  associated  with  the 
experimental  design  (in  this  case,  r<yp2).  In  addition,  the  range  in  the  number  of  facesheet 
plies,  impact  energies,  and  impactor  diameters  considered  in  this  study  increase  the  probability  of 
a  bifurcation  in  failure  modes.  This  may  explain  why  the  response  surface  (equation  9)  did  not 
perform  as  well  in  comparison  to  independent  experimental  data  as  did  the  response  surface 
(equation  7)  developed  in  the  previous  section  for  the  case  where  the  impact  energy  and  impactor 
diameter  were  held  fixed.  It  should  be  noted  that  the  response  surfaces  generally  underpredict 
the  residual  strength  and  overpredict  the  damage  diameter. 
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TABLE  7.  INTERPOLATION  AND  EXTRAPOLATION*  OF  REGRESSION  RESULTS  IN 
THE  SPACE  OF  CODED  MATERIAL  AND  IMPACT  VARIABLES 


X\ 

*4 

*5 

R 

Measured 

Damage 

Diameter, 

D  (in)  r  1 31 

Predicted 

Damage 

Diameter, 

D  (in) 
_ [»3] 

Prediction 

Errors 

D-D 

Measured 
Residual 
Strength, 
Nw  (lbs/in) 

Predicted 

Residual 

Strength, 

Nyy  (lbs/in) 

Prediction 

Errors 

D 

_  .  (%) 

N» 

(%) 

0 

0.03 

-l 

1.000 

1.898 

1.815 

4.4 

1603 

1576 

- 1 _ 

1.7 

0 

0.04 

l 

1.001 

2.046 

2.541 

24.2 

1736 

1498 

13.7 

0 

0.04 

l 

1.001 

2.131 

2.541 

19.3 

1935 

1498 

22.6 

0 

0.05 

l 

1.001 

2.056 

2.551 

23.9 

1790 

1493 

16.6 

0 

0.28 

*1 

1.038 

1.641 

1.802 

9.8 

1894 

1540 

18.7 

0 

0.30 

l 

1.044 

2.851 

2.698 

5.3 

1708 

1424 

16.6 

0 

0.76 

l 

1.256 

2.531 

2.927 

15.7 

N/A 

N/A 

N/A 

0 

-0.98 

l 

1.400 

2.022 

1.739 

14.1 

1985 

1828 

7.9 

0 

-1.01 

-1 

1.421* 

1.603 

1.673* 

4.4’ 

1894 

1766’ 

6.8* 

-1 

0.19 

l 

1.427* 

2.281 

2.130’ 

6.6* 

787 

790* 

0.5* 

-1 

0.22 

1 

1.431* 

2.485 

2.134* 

14.1* 

834 

789* 

5.4* 

-0.73 

1 

1.592’ 

2.283 

1.862* 

18.5’ 

846 

876* 

3.5’ 

-1 

*0.77 

-1 

1.610* 

1.421 

1.908’ 

34.3* 

880 

517* 

41.2* 

1 

0.78 

1 

1.615* 

2.254 

3.867* 

71.7* 

2690 

1562* 

41.9* 

1 

-0.86 

1 

1.655* 

2.291 

2.055* 

10.2* 

2453 

2416’ 

1.5’ 

-0.95 

*1 

1.704* 

1.578 

1.944’ 

23.1* 

964 

514* 

46.6* 

0.97 

l 

1.715* 

2.594 

2.165* 

16.5* 

848 

758* 

10.6’ 

1 

-1.06 

1 

1.767* 

1.724 

1.784’ 

3.4* 

3143 

2529* 

19.6* 

1 

-1.12 

-l 

1.804* 

1.530 

1.483* 

3.2* 

2840 

2796* 

1.6’ 

-1 

-1.17 

i 

1.835* 

1.872 

1.646* 

12.1* 

825 

934* 

13.2* 

0*  =  Extrapolated  values 


2,3  INFLUENCE  OF  FACESHEET  THICKNESS.  IMPACT  ENERGY.  AND  IMPACT 
VELOCITY  ON  THE  IMPACT  DAMAGE  TOLERANCE  OF  SANDWICH  COMPOSITES. 

In  the  final  study  in  this  examination,  the  isolated  effects  of  the  number  of  facesheet  plies  (Y,), 
impact  energy  (X4),  and  impact  velocity  (X6)  on  the  damage  tolerance  characteristics  of  sandwich 
composite  panels  were  investigated  using  the  Box-Behnken  experimental  design.  Here,  different 
impact  velocities  resulting  in  a  particular  value  of  impact  energy  were  obtained  by  varying  the 
impactor  mass  and  drop  height  during  the  impact  testing  (cf.,  [5]).  The  influence  of  dynamic 
effects  on  the  damage  tolerance  properties  of  the  given  sandwich  composites  over  a  range  of 
relatively  low- velocity  impacts  is  being  studied.  Table  8  summarizes  the  low,  midrange,  and 
high  levels  of  the  natural  (Xj,  X4,  X6)  and  coded  (x\,  x4,  x6)  independent  variables  considered  in 
this  effort.  The  core  density  (X2  =  3.0  lb/ft3),  core  thickness  (X3  =  3/4  in),  and  impactor  diameter 
(X5  =  3.0  in)  were  held  fixed. 
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TABLE  8.  NATURAL  VALUES  AND  CORRESPONDING  CODED  LEVELS  OF  THE 
FACESHEET  THICKNESS  AND  IMPACT  VARIABLES 


i 

Independent  Variable 

Natural  Value,  Xt 

Coded  Level,  x, 

Number  of  Facesheet 
Plies 

2  T90/45]] 

-1 

i 

4  r90/45]2 

0 

6  r90/45]3 

+1 

Impact 

Energy 

90.0  in-lbs 

-1 

4 

120.0  in-lbs 

0 

150.0  in-lbs 

+1 

65.21  in/s 

-1 

6 

Impact  Velocity 

96.30  in/s 

0 

127.39  in/s 

+1 

Table  9  summarizes  the  combinations  of  coded  independent  variables  as  well  as  the 
experimentally  measured  CAI  strengths  for  the  15  experiments  used  in  the  regression  analysis. 
Following  the  procedure  outlined  earlier,  statistically  reliable,  second-order  response  surfaces 
were  generated  that  characterize  the  CAI  residual  strength  as  a  continuous  function  of  the  number 
of  facesheet  plies,  impact  energy,  and  impact  velocity.  An  estimate  of  the  CAI  residual  strength 
from  the  regression  analysis  may  be  expressed  either  in  terms  of  coded  or  natural  values  of  the 
independent  variables,  i.e., 

Nyv  =  1820  + 776.9 -x,  -135.65  -x4-  4.930  -x6 

...-97.67 -x?  + 153.9  -x4-  398.9  -x62  (lbs/in)  (10a) 

...  — 119.0-Xj  -x4  -  20.09 -x,  x6  +  21.65  x4  x6 

and 

Nyy  =-1740  +  853.0  -X{  -39.86  -X4  +  77.83  -X6 

...- 24.42  -V,2+  0.1710  -Xl-  0.4127  ■  X]  (lbs/m)  (10b) 

...- 1 .984  -XrX4-  0.323  l-XrX6  +  0.0232  •  X4  ■  X6 

The  constant  term  (bo  -  1820  lbs/in)  in  equation  10a  represents  the  mean  CAI  residual  strength 
for  those  panels  corresponding  to  the  center  point  of  the  design.  Consistent  with  earlier 
observations,  increasing  the  number  of  facesheet  plies,  xi,  from  their  midrange  value  in  equation 
10a  will  result  in  a  linear  increase  in  the  predicted  response,  whereas  increasing  the  impact 
energy,  xa,  and  impact  velocity,  X6,  will  result  in  a  somewhat  smaller  linear  decrease  in  the 
estimated  residual  strength  from  bo.  Note  that  the  quadratic  term  involving  the  impact  velocity 
contributes  significantly  to  the  estimated  response  (i.e.,  the  midrange  velocity  likely  defines  a 
relative  maximum  or  ridge  in  the  predicted  response).  Similar  arguments  can  be  made  when 
interpreting  the  influence  of  the  remaining  quadratic  and  coupling  terms  in  equation  10a  on  the 
estimated  residual  strength. 
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TABLE  9.  COMPARISONS  BETWEEN  PREDICTED  AND  MEASURED  RESIDUAL 

STRENGTH 


Test, 

k 

*i 

x4 

*6 

Measured  Residual 
Strength, 

Nn.  (lbs/in) 

Predicted  Residual 
Strength, 

N„  (lbs/in) 

N>y 

(%) 

31 

+  1 

+1 

0 

2690 

2399 

10.8 

32 

+  1 

-1 

0 

3143 

2908 

7.5 

33 

-1 

+1 

0 

847.8 

1083 

27.8 

34 

-1 

-1 

0 

824.9 

1116 

35.3 

35 

+1 

0 

+1 

1774 

2076 

17.0 

36 

+1 

0 

-1 

1901 

2126 

11.8 

37 

-1 

0 

+1 

786.7 

562.1 

28.6 

38 

-1 

0 

-1 

833.8 

531.7 

36.2 
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2,3. 1  Discussion  of  Facesheet  Thickness  and  Impact  Parameter  Regression  Model  Results. 

Table  9  summarizes  the  measured  and  predicted  CAI  strength  for  each  of  the  sandwich 
composite  panels.  The  difference  between  the  experimentally  measured  CAI  residual  strengths 
and  the  estimated  values  using  equation  10  varied  between  0.2%  and  36.2%  for  the  15  panels 
tested,  with  a  mean  difference  of  13.2%.  The  lack  of  correlation  was  especially  pronounced  for 
xx=  -\  (Xi  =  2  plies).  A  comparison  of  the  data  from  tables  3,  6,  and  9  considered  in  the 
development  of  the  response  surfaces  (equations  7,  9,  and  10),  respectively,  suggests  that  those 
specimens  with  two-ply  facesheet  configurations  often  showed  the  biggest  disparity  between  the 
measured  and  estimated  CAI  residual  strengths.  As  mentioned  previously,  the  compressive 
failure  mechanisms  (facesheet  fracture,  facesheet  dimple  propagation,  etc.)  and  residual  strength 
are  likely  influenced  by  whether  or  not  partial  facesheet  penetration  occurs;  the  nature  and 
severity  of  the  local  damage  also  may  be  drastically  different  for  the  two  cases.  This  is  of 
particular  concern  for  those  impacts  involving  smaller  diameter  impactors,  increasing  impact 
energy,  and/or  thinner  facesheets.  Quadratic  response  surfaces,  such  as  equation  10,  have 
difficulty  predicting  the  CAI  strength  for  those  test  panel  configurations  exhibiting  a  bifurcation 
between  failure  modes  over  the  range  of  input  parameters  considered.  Hence,  the  response 
surface  estimate  (equation  10)  is  likely  inaccurate  for  the  class  of  blunt  object  impacts  considered 
here.  Qualitative  consideration  of  the  response  surface  in  equation  10,  however,  may  prove 
useful  in  discerning  the  influence  of  the  number  facesheet  plies,  impact  energy,  and  impact 
velocity  on  the  CAI  residual  strength.  The  regression  results  may  be  viewed  schematically  in  the 
space  of  coded  independent  variables,  as  shown  in  figure  1 6. 
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FIGURE  16.  RESPONSE  SURFACE  ESTIMATES  OF  THE  RESIDUAL  STRENGTH  AS  A 
FUNCTION  OF  MATERIAL  SYSTEM  AND  IMPACT  PARAMETERS 

Figures  17(a),  17(b),  and  17(c)  show  the  CAI  residual  strength  predicted  using  equation  10  as  a 
function  of  impact  energy  (X4,  x4)  and  impact  velocity  (X^,  xe)  for  the  two-ply,  four-ply,  and  six- 
ply  facesheet  configurations  (or  xi  =  -l,0,  1),  respectively.  Similar  to  the  results  from  the 
previous  sections,  a  comparison  of  the  three  figures  reveals  that  the  predicted  strength  increases 
substantially  as  the  number  of  facesheet  plies  is  increased.  In  addition,  the  maximum  CAI 
strength  for  a  given  combination  of  facesheet  configuration  and  impact  energy  occurs  at  the 
midrange  impact  velocity;  i.e.,  the  CAI  residual  strength  is  somewhat  rate  dependent.  The  three 
contour  plots  are  very  nearly  symmetric  about  a  relative  maximum  or  ridge  in  the  predicted 
response  occurring  in  the  vicinity  of  the  line,  X(,  =  96.30  in/s  (xg  =  0).  The  magnitude  of  the 
estimated  strength  can  decrease  somewhat  as  the  impact  velocity  is  varied  from  the  midrange  test 
value  (x(,  =  0)  to  either  high  or  low  values  (x<s  -  ±1).  Lacy  et  al.  [13]  noted  that  the  predicted 
planar  damage  dimension,  D,  was  also  a  relative  minimum  along  the  line  *6  =  0  for  the  same  set 
of  test  specimens.  Hence,  the  residual  strength  is  somewhat  proportional  to  the  size  of  the  planar 
damage  area.  Sandwich  panel  stiffness  properties,  energy  absorption  capability,  and  support 
boundary  conditions  all  play  a  key  role  in  the  dynamic  impact  response  leading  to  damage 
development  and  loss  of  strength.  This  underscores  the  importance  of  adequately  characterizing 
the  expected  impact  scenarios  when  establishing  a  damage  tolerance  plan  for  sandwich 
composite  aircraft  structures. 
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FIGURE  17.  PREDICTED  RESIDUAL  STRENGTH  (lbs/in):  (a)X,  =  2  plies  [90/45] 
(b)  X\=4  plies  [90/45]2,  AND  (c)  X\=6  plies  [90/45]3 


The  preceding  results  are  in  general  agreement  with  the  experimental  observations.  Figures  18 
through  20  show  the  recorded  CAI  residual  strengths  for  the  test  specimens  corresponding  to  the 
midrange  impact  velocity  (x6  =  0;  X(,  =  2  in),  midrange  impact  energy  (x4  =  0;  X4  =  120  in-lbs), 
and  midrange  facesheet  configurations  (xi  =  0;  X\  =  4  plies),  respectively.  Included  with  each 
observation  is  its  associated  TTU  C-scan  image  and  measured  diameter,  D,  from  reference  13. 
These  figures  show  that  the  experimentally  observed  CAI  strengths  for  a  given  facesheet 
configuration  are  consistent  with  the  regressions  of  figure  16. 
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FIGURE  20.  MEASURED  DAMAGE  DIAMETER  [13]  AND  RESIDUAL  STRENGTH 
FOR  FACESHEET  CONFIGURATION,  Xx  =  4  plies  [90/45]2  (x,  =  0) 


For  the  test  matrix,  summarized  in  table  9,  figure  21  contains  a  comparison  of  the  response 
surface  predictions  of  the  TTU  C-scan  diameter  from  reference  13  to  the  response  surface 
predictions  of  the  normalized  CAI  residual  strength  developed  in  this  examination.  As  noted 
earlier,  the  predicted  normalized  CAI  residual  strength  for  a  given  facesheet  configuration  is 
generally  inversely  proportional  to  the  size  of  the  planar  damage  area,  D.  This  is  consistent  with 
the  observations  of  Tomblin  et  al.  [5]  (cf.,  figure  3).  A  comparison  of  figures  21(d),  21(e),  and 
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21(f)  suggests  that  the  CAI  residual  strength  can  range  roughly  between  50%-70%  of  the  virgin 
panel  strength  for  a  given  facesheet  configuration. 
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2.3.2  Correlation  With  Independent  Experimental  Results. 


Table  10  contains  a  limited  set  of  additional  independent  experimental  results  from  Tomblin  et 
al.  [5],  where  the  combination  of  configuration  and  impact  variables  lie  within  the  spherical 
domain  of  the  test  matrix  design  used  in  this  study.  Response  surface  predictions  of  the  planar 

diameter,  D,  of  the  internal  damage  from  reference  13  and  the  CAI  residual  strength  (equation 
10)  are  compared  to  measured  damage  diameters,  D,  from  reference  13  in  table  10.  The 
measured  and  predicted  CAI  residual  strengths  are  also  compared  in  table  10.  The  difference 
between  the  experimentally  measured  strengths  and  the  interpolated  values  varied  between  1 .4% 
and  23.4%,  with  a  mean  difference  of  13.4%.  Again,  the  difference  between  the  experimentally 
measured  CAI  residual  strengths  and  the  estimated  values  using  equation  10  was  especially 
pronounced  forxi  =  -1  (Xi  =  2  plies)  where  the  residual  strengths  are  greatly  overpredicted. 


TABLE  10.  INTERPOLATION  OF  REGRESSION  RESULTS  IN  THE  SPACE  OF  CODED 

MATERIAL  AND  IMPACT  VARIABLES 
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In  the  final  part  of  this  investigation,  the  coupled  influence  of  the  number  of  facesheet  plies, 
impact  energy,  and  impact  velocity  on  the  CAI  residual  strength  was  evaluated  using  empirically 
based  response  surfaces  for  sandwich  composites  comprised  of  carbon-epoxy  woven  fabric 
facesheets  and  Nomex  honeycomb  cores,  where  the  core  density,  core  thickness,  and  impactor 
diameter  were  held  fixed  in  this  examination.  Lack  of  correlation  between  regression  results  and 
experimental  observations  for  sandwich  panels  with  two-ply  facesheets  suggests  that  quadratic 
response  surfaces,  such  as  equation  10,  may  have  difficulty  predicting  the  CAI  strength  for  those 
test  panels  where  facesheet  penetration  and/or  bifurcation  between  failure  modes  are  occurring. 
Similar  to  the  results  obtained  using  equations  7  and  9,  response  surface  estimates  of  the  CAI 
residual  strength  using  equation  1 0  suggest  that  increasing  the  number  of  facesheet  plies  results 
in  a  significant  improvement  in  the  damage  tolerance  properties.  Regression  results  suggest  that 
damage  development  [13]  and  the  CAI  strength  are  somewhat  sensitive  to  the  velocity  of  the 
impactor.  Midrange  values  of  impact  velocity  resulted  in  strength  estimates  that  were  a  relative 
maximum. 


3.  CONCLUSIONS  AND  RECOMMENDATIONS. 


In  this  investigation,  the  influence  of  sandwich  configuration  and  impact  parameters  on  the 
damage  tolerance  characteristics  of  sandwich  composites  comprised  of  carbon-epoxy  woven 
fabric  facesheets  and  Nomex  honeycomb  cores  was  evaluated  using  empirically  based  response 
surfaces.  A  series  of  carefully  selected  tests  were  used  to  isolate  the  coupled  influence  of  various 
combinations  of  the  number  of  facesheet  plies,  core  density,  core  thickness,  impact  energy, 
impactor  diameter,  and  impact  velocity  on  the  compressive  residual  strength  degradation  in 
sandwich  composites  due  to  impact  normal  to  the  surface  with  relatively  blunt  spherical  steel 
impactors.  The  ranges  of  selected  material  parameters  were  typical  of  those  found  in  common 
aircraft  applications.  Quadratic  response  surface  estimates  of  the  CAI  residual  strength  as  a 
continuous  function  of  material  system  and  impact  parameters  were  developed. 

For  a  fixed  set  of  impact  parameters,  response  surface  estimates  of  the  CAI  residual  strength  are 
highly  sandwich  configuration  dependent.  Regression  results  suggest  that  increasing  the  number 
of  facesheet  plies  results  in  the  greatest  improvement  in  the  damage  tolerance  properties.  In 
addition,  increasing  the  thickness  of  the  core  material  can  produce  a  significant  increase  in  the 
predicted  CAI  strength.  For  sandwich  panels  where  facesheet  penetration  was  not  a  concern, 
residual  strength  estimates  generally  were  inversely  proportional  to  the  size  of  the  planar  damage 
region  obtained  from  a  previous  study.  The  developed  response  surface  for  the  residual  strength 
for  the  above  set  of  impact  parameters  was  sufficiently  accurate  considering  the  scatter  of  the  test 
data  used  to  develop  the  response  surface. 

In  addition,  the  experimental  results  and  regression  analysis  suggest  that  impact  damage 
development  and  loss  of  strength  in  sandwich  composites  is  highly  sensitive  to  the  diameter  of 
the  impactor,  impact  energy,  and  impact  velocity.  An  increase  in  the  impact  energy  and/or 
impactor  diameter  generally  resulted  in  a  decrease  in  the  CAI  residual  strength  for  those  cases 
where  facesheet  penetration  was  not  likely.  Moreover,  response  surface  results  indicate  that  the 
damage  formation  and  loss  of  strength  is  somewhat  sensitive  to  the  velocity  of  the  impactor. 
Sandwich  panel  stiffness  properties,  energy  absorption  capability,  and  support  boundary 
conditions  all  play  a  key  role  in  the  dynamic  impact  response  leading  to  damage  development. 

For  a  given  range  of  sandwich  configuration  and  impact  parameters,  response  surfaces  may  be 
used  to  estimate  the  size  of  the  ensuing  impact  damage  as  well  as  the  degradation  in  residual 
strength.  Special  care  should  be  taken  when  estimating  the  CAI  residual  strength  for  cases  where 
the  range  of  independent  variables  spans  a  potential  bifurcation  between  failure  modes.  This 
occurs  primarily  for  thin  facesheets  and  sharp  impactors.  A  response  surface  for  two  diameters 
and  two  facesheet  thicknesses  may  be  too  narrow  in  scope  to  be  of  any  practical  use.  In  general, 
the  response  surfaces  for  the  compressive  residual  strength  developed  here  can  be  of  qualitative 
use,  but  are  not  as  useful  as  those  developed  for  damage  resistance. 

To  fully  develop  the  concepts  outlined  here,  a  detailed  risk  assessment  establishing  impact 
scenarios  that  represent  a  viable  threat  to  principle  structural  elements  is  required  prior  to 
establishing  a  damage  tolerance  plan.  Such  an  assessment  would  likely  involve  specifying  the 
upper  bounds  on  the  impact  parameters  (impact  energy,  impact  velocity,  and  impactor  diameter) 
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that  represent  a  realistic  threat  to  a  given  structural  component.  It  is  important  that  the  structural 

integrity  is  not  seriously  degraded  by  any  impact  damage  falling  below  the  threshold  of 

detectability  using  standard  detection  techniques. 
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